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I. INTRODUCTION 


Aa  FURFOSE 

The development of digital transmission techniques in 
military communications systems has generated a renewed 
interest in tropospheric scatter. Digital troposcatter 
radio systems have required changes in the present 
planning and engineering design methods. The military 
design engineer has the responsibility to predict system 
performance and reliability within limited planning 
time. Under certain tactical situations, the urgency to 
restore vital communications has iore. the engineer tea 
provide  zeveral alternative designs f OF а Single 
communications link. These requirements have justified 
the need for an efficient, computer-aided design method 
to accomplish the task. The purpose of this thesis 
research 1= τω develop a tropospheric scatter 
communications system design program. The objective is 
to design а microcomputer-based program which will 
predict the path reliability and system perfarmance far 


troposcatter  radiolinks. This program will introduce 


several improvements Over the manually generated, 
nomograph-based design methods. Specifically, the 
program will estimate the height gain at the average 
system antenna height for specific elevated tropospheric 
ducts. Other features include a terrain profile pict, 
radiosonde weather data analysis, and the calculation of 


the probabilty of error for binary and orthogonal digitai 


Signaling techniques. 


BACKGROUND 
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Tropospheric scatter communications systems pords 
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microwave energy beyond line-of-sight (EOS) ter. Paver th 


It 


horizon" һу taking advantage of the refraction and 
reflection phenomena in a section of the earth's 


atmosphere called the troposphere. Typical military 
troposcatter systems use transmitter power outputs from i 


t 10 kW, parabolic type antennas, and sensitive 


Ü 


broadband FM receivers with front-end noise figures (MF? 


г 


jgtween 2.0 — 43.0 dE. These systems are mas orten 
employed in military applications because of the limited 
available channel capacity CRef. lipp. 252-2541]. Severe! 


system advantages are summarized as follows: 


і. Fewer terminal and/or relay stations are required 
to cover the same transmission path as compared 
to tactical microwave line-of-sight radia es 


2. Reliable multichannel communications can be 
installed across long distances of hostile or 
inaccessible terrain. 

5. Standard transmission ranges are suitable for the 

d 
tactical military field environmemt for radiolinksz 
from 20 to 200 miles. 

A fundamental design parameter for troposcatter links 
is the allocation of allowable channel degradation toa 
each link making up the total communications system. in 
analog scatter systems, the quality parameters are system 
availability and signal-to-noise ratio in the σι... 


channel. However with the deployment at digit = 


troposcatter radio systems the perrorma n -TEs the 


į 


digitized voice channels under fading conditions is 
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bel 


- 
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fu 
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Significantly different than the performance of og 
yoice systems under the same conditions. The primary 
measure of transmission quality for a digital system is 


its error performance. There are two main sources TT 


transmission errors: (a) long-term error rate which will 


(I 
usd 


Occur because Of equipment degradation, chann 


Cr 


interference, and long-term power fading; and i 


multipath тай“ та Chet. οσο uu 


C. RELATED WORK 
1. Performance Models for Troposcatter Links 

Monsen CRef. 31 has derived performance models 
for analog  FDM/FM and digital quadrature phase shift 
keying (QF'SE) systems on troposcatter communications 
links. The analysis included the etfects of signal level 
variations and multipath delay dispersion. His 
performance criterion was the outage probability rather 
than either average error rate or median signal-to-noise 
ratio (SNR) values. Outage probability performance in a 
digital system was derived for two different modems: thea 
Decizion-Feedback Equalizer (DFE) and a transmitter 
time-gating technique, the Distortian Adaptive Receiver 
(DAR). Analog system performance was determined for tha 
percentage ' of the time that the voice channel SNR 
including thermal noise and multipath delay dispersion 
effects were below a specified threshold. 

The performance test results for a Distartian 
Adaptive Receiver (DAR) was presented in a paper by 
Zawislan CRef. 4]. The critical problem with digital 
troposcatter is intersymbol distortion produced ἂν tha 
multipath dispersion of the fading channel. Та resolve 


this problem, the distortion adaptive receiver Was 


implemented. The DAR modem Та OFS modulation with 
adaptive matched filter demodulation. A complete 
functional description was explained in the reference. 
The advantage to this approach was that it did not 
require equalization at the receiver to correct 
intersymbol interference and it provided near optimum 
performance using the adaptive matched filter receiver. 
Typical bit error rate performance for the DAR 
on a troposcatter channel is shown in Figure i-i For 
different values of rms multipath delay dispersion. Fer 
low multipath spread, the fading on each diversity 
channel has a Rayleigh distribution and the pertormance 
is indicated by the "flat fading" marked curve. Mote 


that the performance improves rapidly with increasing 


с 
ч 


multipath spread and, at a іо error rate a іс а 


ri 


improvement is achieved for multipath curve CBI]. However, 
as the multipath spread increases (with increased pati 


length), а irreducible error rate occurs as hawn in 


Ifi 


curve [C] due to the transmitter time-gate limitation. 
To compensate for this problem, the DAR is moditied Бу 
utilizing a dual frequency pulse wavefarm. This 
technique 15 also used in the modem employed By the 


AN/TRC-170 tactical digital troposcatter radio set. 
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Figure 1-1 Distortion Adaptive Receiver Ferformance 


(After Zawislan, Ref. 4:p. 3) 


т» Tropospheric Elevated Duct Models 


several mathematical models have been developed 

to calculate the ducted signal field strength based on 
the refractivity profile and by modeling the tropospheric 
duct as a laterally homogeneous-layer atmospheric 
waveguide [Net о рр οσους à detailed mathematical 
analysis was presented by Marcus and Stuart CRef. ó]. 
These authors have developed an effective modal solution 
computer program, called "DUCT", which requires a large 
processing capabilty of a main-frame computer CRef. 7]. А 
version of this program, called "FDUCT", is available at 
the Naval Fostgraduate School Computer Center (CRetf. δη 
апа LASt. S600. cad 


Weston CRef. S:pp. 38-701] has proposed a 


statistical model, based on  "PDUCT" predictions for 
selected ducts of interest. The height gain data, for 
these ducts, was used to derive model coefficients. 


These unique coefficients were stored according to their 
respective transmitter and receiver heights in a matrix 
form. Thus, given a receiver height and range from the 
transmitter for a selected duct/frequencv data point, the 
coefficients for each receiver height are obtained апа 


used ta produce the corresponding height gain curve. 


This approach required a significant number of main-frame 
computer "FDUCT" computations to establish height gain 
data for all the ducts of interest. 

In December 1983, a microcomputer-based Bro- 
gram, | called  "MIMIDUCT" Version 1.1, based оп Enorr’s 
mathematical model was developed by Nagel, CRef. 91, at 
the Naval Fostgraduate School to calculate elevated 
ducted signal levels at various frequencies. This 
program used either historical radiosonde data or current 
elevated duct information [Ref. 10]. The radiosonde data 
represented the atmospheric temperature, pressure, sand 


humidity at specific altitudes over a five year recordin 


Ц) 


period. This model was valid for the case where tie 
transmitter is at the optimum duct coupling height and 
the receiver is located either below, within ar above the 
duct. The optimum coupling height is the altitude above 


the surface where the gradient Of the modified refractiv 


Ш 


index becomes negitive. 


II. TROFOSHERIC FERFORMONCE CONS DE 


a. SCATTER FROFAGATION 
te General 


At frequencies above 30 MHz three propagatian 


mechanisms can carry energy beyond the horizon: (a) 
variations in the refractive index in the troposphere can 
scatter radio energy, (b) Ν᾽; abruet 
changes in the refractive index can cause reflection, ana 
(C) atmospheric regions of negitive modified refractive 
index gradients can introduce ducting. Refractive index 


and tropospheric ducting will be thoroughly discussed in 


the following section. The forward scattering aft radia 


mechan 


[54 
П 
a 
ib 


= 
С 
= 
Πι 


signals 1s the most dominant propagat 
frequencies of 0.5 to 10 GHz. Lmerc lias 
РЕ Received Scattered Field 

The index Gf refraction denends on pressure, 


humidity, and temperature. Slight variations in these 


quanititi 


{1 
iñ 


caused hy atmospheric turbulence, wap bd 
produce random flucuatians i the refractive index. Wen 


an electromagnetic Wave propagates through tnis 


inhomogeneous medium, energy will be scattered out from 
the original incident direction. The 
turbulent-scattering theory, [Ref. 12:p. 3451, has shown 
to a first approximation that the index of refraction 
fluctuations can be replaced by a model of so-called 
Blobs", inhomogeneities of different dielectric 
constants randomly distributed. If these "blobs" are in 
the common volume formed by the transmitter and receiver 
antenna beams, the complex received field сап be 


described by [Ref. 13:p. 14051: ` 


f m : 
207 Σ jo. 
Re = DE e 1 
=] 


where m is the number of "blobs" in the scattering 
volume, Ar is the amplitude, and ф is the phase of a 
wave MS red by a single "Blob". Assume m to be very 
large and the blobs are spherical апа uniformly 
distributed through the scattering volume. Then the 
phase difference of the waves scattered by the 


inhomogeneities at the top and bottom of the volume will 


Def Ref. 15:р. 147]: 


"PEDEM 


| 
۲ 
t 


17 


where y transmitter, receiver antenna elevation angle 


h = mean thickness of the scatter volume (т) 


À 


transmitter wavelength (т) 
The variables in Equation 2-2 are described in Figure 2-1. 
By observing the volume geometry, it can be determined 
that more "blobs" exist in the central volume region than 
at the top or bottom. [his indicates a  non-uniform 
"Blgcb' gristributisnmn. 

Many phase variations will occur over many 


oh 


ni 


se cycles of length 27 . Тһе antenna elevation analis 


Will be less than £ degrees for most system 


ul 
HT 
" 
ek 
rr 
ср 
i 


ratio of scatter volume thickness to operating wavelength 
will σε very large. These conditions make the scattered 
phases Ey Yee ml distributed from zero ta m 
Кю» arpe 1481. 

The amplitude distribution will now be deter- 
mined. The real and imaginary components of the complex 
received field, Equation  Z-1, were resolved into two 


random variables X and Y. The number of blobs m, 1 


ill 


assumed to be large, so by the Central Limit Theorem 
Foth X and Y will approach a normal “gaussian” 


πο ει Επιστ Finally it was shown that these variabies 
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Figure 2-1 Scatter Volume Geometry 


(After Beckmann, Kef. 11:p. 147) 


are also uncorrelated апа independent. Then their 
two-dimensional normal distribution can be transformed 
into the Rayleigh distribution. 
Turbulence, mixing, and wind continuously 
change the positions and structure of the inhomogeneities 
making the random terms of Equation 2-1 functions of 
time. The reality of nature is the inhomogeneities do 
not vary isotropically because the index of refraction 
changes more rapidly with altitude than with distance. 
There can also be a stratification or layering of 
inhomogeneities near a particular height (or heights) 
à 

that changes the assumed spherical shape of the so-called 
"blobs". Under these conditions the uniform phase 


distribution may change to an unknown distribution. 


E.  TROFOSFHERIC DUCTANG: 

The index of refraction, n, of air is defined as 
the ratio of the velocity in a vacuum of electromagnetic 
(EM) radiation to the velocity in the medium. A 
convenient parameter is refractivity, N, defined as 


Ке. ел λα, 


1 
~ 


Another parameter called modified refractıvity, M, is 


defined as (Ref. 10:p.91: 


Me = Πο nh 


where h is the altitude in meters above the surface. The 
modified refractivity accounts for the curvature of the 
earth so the presence of ducting can be easily determined 
by observing the M-gradient on the M versus height plot. 
Refraction of incident radio waves across a 
discontinuity of refractivity is described by the 
principles of Snell's Law. It is important to remember 
that the wave "bends" towards the higher value of 
refractivity, and the more dense a material the higher 
its n. Since the density of the atmosphere decreases 


with height, we expect that a wave will bend back 


downward from a geometric straight path. Whenever the 


refractive index decreases sharply with height, 
radiowaves can be trapped апа experience low-loss 
propagation for long distances. This condition is known 


as tropospheric ducting [Ref. lisp. 29]. 
The following conditions must be satisfied for a duct 
ο σεςειιγ- (1) the modified refractive index gradient 


shall be equal to or more negitive than O M-units/km, and 


rJ 
NS 


(2) this cradient should continue over a height of many 
wavelengths. The important duct parameters are the duct 
thickness, D, the intensity, M, and the optimum coupling 
herght. BH . A piecewise linear approximation ta the 
opt 
modified refractivity (tri-linear) profile for several 
types of ducts are shown in Figure 2-2. There are three 
types of ducts: cl; surface or ground-based ducts, also 
called evaporation ducts when formed over water, Figure 
с-га (2) surface-based ducts from elevated refractive 


layers, Figure 2-2b; and (3) elevated ducts from elevated 


refractive layers, Figure/’s) 2-29с and d. Noto πας -- 


m == 


positive M-gradients are assumed at 118 M-unitzs/kim which 
corresponds to a standard atmosphere. Once the slopes 
are identified, the important duct parameters are quickly 
determined thet. етра οτι 

Tropospheric ducts more often occur as qround-based 
ducts because of both, evaporation and advection. 


Evaporation of water vapor from the surface af the sea 
may create a zone of high humidity (i.e. high retractive 


index? below a region af drier air. Advectian, derinea 


movement of ane air type over another, may cause 


μ.ο 
c 
{з 
it 
ci 


πως dry air (from the land) to be blown over co 


air, producing a region of low refractive index above 


у 
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region of high refractive index. Such a duct παν also 
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Figure 2-2  Tropospheric Duct Descriptions 


(After Marcus, Ref. 7:p. 2-9) 


I. 
L. 


form when warm dry air is blown over cold qraund. 
Radiation cooling can produce temperature gradients which 
cause ground-based ducts. Air next to the qround becomes 
ceoler and the duct becames thicker as the night 
continues. 

When morning solar heating warms the air next to the 
Ground, a region of rapid decrease of refractive index 
with height produces an elevated duct. However these 


ducts quickly disappear because continued ground heating 


increases convection mixing and destroys the stable 
elevated layer. Elevated ducts may farm far several dave 
by 8 subsidence inversion. Hot air rises at the center 
Gf a high pressure region and spreads out horizontally. 


cooling as it slowly descends. This produces a boundary 


-- 


with the slightly colder air near the surface. he 
increasing temperature with height at the boundary formsa 


the subsidence (temperature) inversion. Changes in 


temperature and/or humidity may cause related changes in 


- 


the refractive index within the houndary interval 
rs ЕЕЕ е а 
Figures  Z—i2a thru d illustrate propagating ja 


سر 


within ground-based ducts. Rays leaving the transmitt 
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Figure 2-3 Fropagation Effects of Ducts 


(After Hall, Ref. 9:p. 31) 


2 


either travel upward or downward. If the M-gradient 15 
negitive and assumed uniform within the duct, rays 
beginning at elevation angles below the critical angle, 


defined as (CRef. lisp. 311: 


0 = ү AM X QN (2-5) 


C 


will strike the surface, and those above the critical 
angle will leave the duct. If ground reflection is 
neglected, Figure 2-3a, a shadow region depicts minimum 
propagated energy. In the usos where surface 
reflectivity is high, Figure 2-Sb, the rays launched 


within the critical angle range will travel beyond the 


radio horizon. Figure 2-3c decribes the situation when 


the M-gradient is not constant within the duct. If Be 

transmitter is positioned in the height region h , then 
i 6 я 

rays within the t 0 | = = e AM, IOS range will 


remain in the duct with bounces as depicted in Figure 


Πα. Similarly rays transmitted within 





зеке. == (2 AMS A O ee E are trapped within the 


height range h . Finally а surface-based duct from an 
b 
elevated layer over rough terrain is illustrated in 


Figure 2-3d. In this case the refractive "bending" takes 


place at the top of the duct [Ref. Ἱτ πο LE 


For elevated ducts formed by advection or subsidence, 
the position of the transmitter and receiver relative to 
the optimum coupling height, H quill influence the 

opt 

propagation effects. The means in which energy enters ar 
leaves an elevated duct can be described by the duct 
acting as a "leaky" wavequide. Energy is "leaked" ar 
coupled into the duct from the transmitter, and "leaks" 
aut а= the energy propagates alang the ΕΤ 
NULL 11:р. 26). 

Бесац=е of the non-permanent characteristic = 


elevated ducts, their effects are seldom an influence to 


troposcatter links, especially if they form above the 
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ШЫ с Г ише SUC ERE presence of troposph 
ducts can degrade the overall performance of troposcatter 


systems бу changing the predicted transmission loss. TH 


10 


term that will change the total path loss iz the duct’ 


μὴ 


height gain, which i 


ui 


derived in Chapter IV. 


See MULTIFATH CONSIDERATIONS 


The multipath fading model for a tropospheric scatter 


channel produces received signal ‘fading. The received 
ΓΙ 15: consists of the sum of a large number Gf 
time-variant, complex vectors having amplitudes anda 
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caused by randamnly ime-vaeriant 


phases variations. At times the received signal vectors 
add destructively to decrease the mean received signal 
amplitude. While at other times, the vectors add 
constructively, so the received signal is large. Thus 
the amplitude variations, or signal fading, are due to 
the multipath characteristics of the tropospheric 
channel. The channel can be modeled as a zero mean, 
complex-valued gaussian process, with the envelope of the 
instantaneous signal level being Fayleigh-distributed. 
This Rayleigh fading channel describes the short-term 
fading. When there аге fixed constant regions ar 
refractivity or stratified refractive layers in the 
vicinity of the cammon scatter volume the channel cannot 


De modeled az having a zero mean. In this case, the 


Rayleigh distribution does πος apply and the chanmol 


ui 


aporoache 


a Rice distribution (Ref. it:pp. 456-4581. All 
performance predictions in this study will assume 
Rayleigh "short-term" fading as the channel model. 

Channel performance will be degraded during periads 


of severe multipath fading. With digital svstem 


Iñ 


/Dice user is unaware of any increase in hackgrcound noise 


ii 


until the FCH (Pulse Code Modulation? outage threshaoid i 


broken. Once this threshold is passed the complet 


it 


iti 


rt 


- 
iD 
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αι P i cia s) sssi i will be unusable cude torrak e] 


characteristics of fade Outages for a typical digital 
troposcatter circuit can experience three (5) primary 
categories of fade outage. The first category occurs 
when the voice user is subjected to a single fade outage 
of duration less than 200 milliseconds. This outage will 
be hardly noticed. The second category of outage will 
have an outage duration ranging from 1/5 second to ° 
seconds. The user will detect this distortion but will 
continue to communicate following the outage. When a 
recurrence of short duration outages take place (e.g. 2 
to 4 outages per minute) annoyance rather than total 
disruption will occur. The third category are fade 
Outages that exceed a subjective level of user patience. 
The Defense Communications Agency (DCA) [Ref. 151] has 
epecified five (5) ranges of fade outage conditions, 
refer to Table I . DCA has defined the fade outage in 
terms of a diversity signal-to-noise ratio threshold 
corresponding to a E bit error probability. 

Techniques used to counter multipath propagation are 
frequency | diversity, space diversitv, amplitude 
equalization, and channel equalization. If a modulated 
signal is simultaneously transmitted over the same 
troposcatter radiolink on two or more frequencies, the 


correlation between the individual received signals will 


Qutage Range 


II 
III 
IV 


NOTE: 


Range 


11] 
IV 


TABLE I 


VOICE PERFORMANCE CHARACTERISTICS 


FOR 


DCS  TROFOSCATTER LINKS 


Criteria Outage Probability 
See Note See Note 
_4 
O.2 sec.< Outage <3 Sec. 7.5105 
S sec.< Outage <1 min.  7.5X10-3 
= < Outages/min. < 5 2.5X10-4 
о < Outages/min. 1.0X10 


Voice Performance Description 


Qutages with adequate fade margin. 
Outages with adequate fade margin 
and high frequency of occurrence. 
Call disruption possible. 

Marginal fade margin. 
Unavailability of circuit. 


be small. This method of siqnal diversity is called 
frequency diversity. The important advantaqe cut 
frequency diversity 15 that it requires a single antenna 
at each site. But the need far additional frequencies 
can increase the probability of ca-channel interference 
among other operating transmitters. 


Uncorrelated short-term fading can also be achieved 


Ji 
>. 
й 


Dv separating the receiving antennas in space. Thi 
known as space diversity. Horizontal and vertical 


polarization can “be used to distinguish between twa 


fwd 


epace-senarated signals. However horizontal and vertica 


i 
1 


polarization da not provide a satisfactory degree 


[J 


nancorrelation of signal fading far efficient diversity. 


The multiplicity of signals provided DY these 
diwversitv methods must he cambined. The 
diversity-combining techniques are classitied inta: #@ 
selection or switching; (2) cambining CeCe Ss е ae 


weighted combination of received available signals: and 


(. 3) a cambinatian of selection and combining. 
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rr 
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election techniques the diversity channel 


are scanned until ane is found whose level exceeds a 
selected threshold. This may nat necessarily select the 
best available signal. In the cambining techniques. all 


diversity channels аге simultaneously monitored and 
equally weighted. This is called  equal-gain combining. 
In maximal-ratio combining, the weighting factor of each 
channel is automatically adjusted to yeild the maximum 
signal-to-noise ratio for the total of all the diversity 


channels CRef. 12:pp. 453-454]. 


Amplitude equalizers are designed to proper: 
equalize the propagation channel far minimum phase 
fading. Channel equalizers balance the channel for 
amplitude. and mutlipath delay  distrotion. They are 


typically adaptive transversal equalizers that consist of 
tapped delay lines (TDL) with  tap-weight multipliers 
fi.e., amplifiers or attenuators); and control circuitry 
that adaptively vary the tap-weights in response to 


temporal channel variations πε, δρ. DI DENM 


ΕΣ 


III.  TROFOSCATTER SYSTEM DESIGN PROGRAM 


A. GENERAL FROGRAM ALGORITHM 
The main program, named "TROPO", 15 presented as 

Appendix B. Figure 5-1 illustrates the program flow and 
the primary computational program modules. The program 
development and mathematical approach for the Radiosonde 
Data Analysis and Height Gain modules are described in 
Chapter IV. The remaining modules are formulated in this 
chapter. A program tutorial and compilation instructions 


are contained within the program. 


F. PREDICTION QF FATH LOSS а 
D. General 
[he basic median transmission loss will be the 


sum of several additive losses, expressed in decibels, 


ο ες 171: 


= s Í К+ - = vlc (5-1) 
Lip ES E Lo a ET τ ο. 
where L = free-space propagation/scatter loss (dB) 
s 
L = kEnife-edge diffraction loss (dB) 
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Module Block Diagram 


"TROFO" 


Frogram 
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TA 
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Figure 
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L = medium-to-aperture coupling loss (dB) 


L = atmospheric absorption loss (dB) 
Е, = waveguide/connector loss (dB) 

5 = transmit antenna gain (dB) 

E = receive antenna gain (dB) 

но = height qain (dB) 


E. surface Refractivity 


An adjustment to the average surface refrac- 


tivity N , refer to Figure 3-2, is made for the elevation 
o 
at each terminal site. The adjusted surface 


metmeactivity, N , is (Ref. 18:p. 2-121: 


5 


= i =U a 3222D (Ж— 2) 
N. N exp(-0.0 ο) 


where N minimum monthly average refractivity 
О 
h 


= 


average antenna height (kft) 


If the surface refractivity at each site is significantly 
different, an option to calculate the respective N for. 
each site can be selected and an average path E мс 
refractivity can be calculated. The average antenna 


height 45 calculated by averaging the transmit and 


receive antenna heights. 
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Monthly Surface Refractivity 


(After Panter, Ref. 12:0. 375) 


(Sea Level) 


3. Effective Earth's Radius 
Changes in the surface refractivity will effect 
the radio-ray path curvature as the wave propagates over 
the earth. To accurately represent the ray-path an 
effective earth's radius 1 55 calculated as 


BEef. I2:p. 5741: 


а = a. Ë - 0.04665 exp (0-005577) s (3-3) 
where a = actual earth’s radius (6370 km) 
О 
4. Terrain Frofile Flot 
Several methods are available to plot a 
troposcatter system terrain profile. The most 


fundamental method is to plot the successive path terrain 
elevations along the great circle path. opecial 4/3 
earth plotting graph paper is required for this method. 
Alternatively, computer graphics technigues which obtain 
terrain information from topographical databases, can 
rapidly plot the profile. 
The program provides a dot-matrix printer plot. 

The terrain profile is linearly plotted by modifying the 


terrain elevations, in meters, to include the effect of 


the average curvature of the radio-ray path and the 
earth’s surfave. Elevations, D o, of the terrain are 
manually obtained from Conor eph maps and tabulated 
versus distances, x , from a selected reference 
terminal. The T data points are keyboard entered 


into the program. The modified elevations are computed 


as [Ref. 12:p.380]1: 


2 
х і 
ү. = h, _ (3-4) 
Б 2а 
where a = effective earth's radius (Equation 3-3) 


Figure 3-3 illustrates a typical linear terrain plot. 


= 


o9. Calculation of the Anqular Distance 
The terrain profile can now determine various 

path geometries. The three (3) path configurations 
considered are: 

a. Smooth Earth Horizons at Both Terminals 

b. Obstacle Horizons at Each Terminal 

с. Smooth Earth and Obstacle Horizons 
The terrain geometries are depicted in Figure(s) S3-4a, 
5-46, апа 3-4c. The respective terminal take-off angles 
are calculated τῶι the predicted path type as 


(het. 121855. 55515 71: 
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Figure 3-3 Linear Terrain Frofile Plot Example 


(After Panter, Ref. 12:p. 382) 


41 





Figure 3-4a Smooth Earth Path 


(After Fanter, Ref. 12:p. 585) 





Figure 3-4b Near Obstacle Path 


(After Panter, Ref. 12:p. 386) 





Figure 3-4c Near Obstacle/Smooth Earth Fath 


(After Fanter, Ref. 12:p. 585) 


44 


a. Type I Radio Horizon/Kadio Horizon 


2ah,. 
o е ———— ευ οσα) 
а 
2ah 
Ө n rs (Sab) 
er 
a 
where h , h = transmitter, receiver terminal elevation 


ts rs 


b. Type II Obstacle Horizon/Obstacle Horizon 








_ Ес Tie 
Pet 7 а ST. (3-5c) 
UE 
h = 
С ie hrs dir 
Өег Е F (3-Sd) 
di 2a 
r 
where h zd = transmitter obstacle elevation,distance 
itg lt 
h ΤΩ = receiver obstacle elevation, distance 
irg lr 


C. Type III Obstacle Horizon/Radio Horizon 


Ө = (Same as Equation 3-Sc) 
et 
0 
er 
The path type combination may be reversed to satisfy a 


(Same as Equation 3-Sb) 


Radio Horizon/Obstacle Horizon configuration. Referring 


to Figure 3-53, the angles a апа fj are calculated as: 
оо оо 





αλ ο ο νο, 


Figure 3-S Typical Fath Geometry 


(After Ref. 18:p.8-6) 
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E = ies eee + Oot (5-58) 
OO 2a d 
h - h 
Е rs ts (3-5f) 
Boo š > р бег 
2а а 


These angles are modified by correction factors, Aa, and 
AB, to allow for the effects of а non-linear 
refractivity gradient  [RHef. 12:p. 3831]. The correction 
factors can be obtained from Appendix A, Figure A-7, 
however for most transhorizon, “over the horizon", paths 
these factors are negligible. 

The angular distance (often called scatter 


angle) is: 


INCUN Beet 8 PS (5-5д) 


The ratio Gago and Boog defines the path symmetry factor 


Beet. 18:p. 4-71: 


Cos 
D 


The following equation will calculate the heigth of the 
intersection point of the transmit and receive antenna 
Deams. This result will estimate the bottom of the 


common scatter volume (Refer to Figure SSO) as 


47 





Figure 3-6 Common Scatter Volume Dimension 


(After Du Castel, Ref. 19:p. 145) 


48 


Eef. 18:p. 8—81: 


s d8 
h = ην 
where S = symmetry factor 
Du Castel, [Ref. 19:p 146], approximates 
the dimensions of а typical common scatter volume. The 
volume height is: 
Оо 
H> 2h А 2 = 0.3h Coan) 
О 0 О 
О 


where 6 = 1.1560, 


> scatter angle 


4 
The maximum longitudinal dimension, L, 15: 


= 0-150 


where d - great circle path distance (т) 


of the scatter volume is 


The center of gravity 


approximated as: 


H Sie ο Stl 
cg О 


Ge Diffraction Loss 
Fropagation paths having a common obstacle 
horizon, such as a mountain ridge, can be referred to as 
an obstacle gain path. It 15 assumed that the obstacle 
will introduce additional path attenuation. However, the 
angular distance may be reduced because of the changed 
path geometry from the obstacle. The possible reduction 
in the scatter loss may be offset by the increased loss 
due to diffraction over the obstacle. In some situations 
the common obstacle may be visible to both terminals, and 
the path loss might be less than the smooth earth path 
loss. The International Radio Consultative Committee 
(C.C. ο, has developed the following formula for 


diffraction loss relative to free-space [Ref. 20:р. 1701]: 


ta ta 3 
| tan 8, tang, 





where d = transmitter to obstacle distance (m) 
a 

d 

b 


When the take-off angles are less than 10 degrees and 


receiver to obstacle distance (m) 


d is greater than 2b then Equation 3-5a can be 
b а 
approximated by: 


La 2 201og,), | 276 “x (Ξ- 68) 


7.  MWorst-Hour Loss 

Seasonal annual-to-worst month path loss 
variations can be determined from a knowledge of the 
annual changes in surface refractivity of the atmosphere 
over the path. The C.C.1.R. recommends a loss variation 
of O.2 dB (U.S.) and O.5 dB (Europe) per unit change of 
refractive index. 

The worst-hour median loss can be derived by 
assuming a log-normal distribution during the month. ГЕ 
was shown by Fanter, [Ref. l2:p. 4011, that on a 
log-normal distribution, the 99.9 percent point can be 
approximated by the value 3.10 in decibels below the 
median, where c is the standard deviation of the 
log-normal distribution. The worst-hour median loss can 
be determined as: 

Worst-Hour Median Loss = Median Annual Fath Loss 

* Difference of Annual-to-Worst-Month Median Loss 

+ $3.1 QC wm 


where O = standard deviation of the worst-month 
wm 
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Figure 2-7 Annual and Winter Monthly Standard Deviation 


(After Fanter, Ref. 1Z:p. 359) 


Ch 
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distribution obtained from Fiqure 35-7. 


e Aperture-to-Medium Coupling Loss 


The parabolic reflector microwave antenna 
gain equation, [Ref. 12:р. 1051, can be expressed in 


decibels as: 


G = 20108, f + 20103, و‎ 52.6 OB (° 7) 


where D = aperture diameter (feet) 
An illunimation factor of 0.34 was assumed to derive 
Equation 3-7. It would appear that Equation 3-7 depicts 
an ever increasing power gain as the antenna aperture 
area increases. However the power received by an antenna 
does not increase linearly with an increase in antenna 
diameter, D. This effect is called aperture-to-medium 
coupling loss or loss in antenna gain [Kef. 12:p. 3621. 
Aperture coupling loss in troposcatter systems 
is caused by a non-planar wavefront due to atmospheric 
irregularities, and a geometric effect due to the 
decrease in the scattering properties with height inside 
the scatter volume. An incoming wavefront consists of 
many plane waves, each arriving at a different angle from 
the scatter volume. If the arrival angle range variation 


but 
ted eat 


is much smaller than the antenna beamwidth the wavefront 
will appear nearly plane. If the common volume is much 
wider than the receiving antenna s beamwidth a wider 
angle range will result, and the wavefront will appear 
non-planar. 

Basically, coupling loss can be explained as 
limited antenna pickup, as compared to the effective 
scatter volume dimensions and the antenna 3 dE 
beamwidth. Figure 3-8 compares aperture coupling loss 
results between several authors. A unique constant curve 
is presented by the C.C.I.R. LRef. 21:p. 1451. Ths cu 


is independent of the scatter angle and is written as: 


Lo = 0.07exp | 0.055 (Gy + 6.) | (S 7 


where G ,G = transmit, receive antenna gain (dB) 
T R 


This empirical formula gives a high coupling loss and 
will not be used in the program. The proposed empirical 
Curve appears as the average of several different 
formulas and will be used as a conservative estimate for 


the aperture coupling loss. 
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Figure 3-8  Aperture-to-Medium Caupling Loss 


(After Levin, Ref. 21) 


9. Combined Free-Space/Scatter Loss 


Yeh, tKef. 221, has derived the following 
formula to calculate the combined free-space and scatter 


loss, in decibels: 


= Έτ = (31709 
L. 30103 م‎ H 20 log; 94 ΕΠΕ 0.211. 310) + 597 
where f = operating frequency (MHz) 
d = great circle path distance (miles) 
Ө = Scatter Mole (degi ees) 
N = surface refractivity 
s 


10. Waveguide/Connector Loss 


The waveguide attenuation factor was derived 
for standard rigid waveguide. At an operating frequen 
of 4.5 GHz the waveguide loss will be approximately г 
dB per 100 feet  LhNef. 18:p: 7/7- rF: Each waveguide 
connection will introduce an additional 0.06 dH per joint 
Екей 1:5. 2101- 

li. Absorption Loss 

Rainfall, snowfall, and fog produce atmos- 

pheric absorption loss which depends on the amount of 


moisture and on the frequency. Figure 3-9 was used to 
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Figure 3-9 Estimated Atmospheric Absorption 


(After Ref. 18:p. 2-43) 


estimate the rainfall attenuation on a particular 


transmission path at four (4) rainfall rates 


[Ref. 18:p.2-—431. 


C. DESIGN FARAMETERS 
і. Carrier-to-Noise Ratio 
The performance of troposcatter communica- 

tions circuits are determined by the minimum acceptable 
ratio of hourly median carrier signal to thermal noise 
for a type of modulated signal. This ratio, expressed in 
decibels, is called the carrier-to-noise ratio (CNR). The 
following equation is used to determine the received 


carrier power level: 


P, (dBw) - Ly, (dB) 


Р, (ави) 5-10) 


where Г transmitter power output (dEU) 
t 
Ё 

1 


The receiver noise threshold level is written as: 


total median transmission loss (Equation 3-1) 


p (dBW) = -204(dBW) + NF(dB) + l0log), By, (3-11) 


n 


where —204 dEW - thermal noise constant 


NF receiver noise figure (dB) 


B receiver IF bandwidth (Hz) 
IF 


Finally the carrier-to-noise ratio, [Ref. 10:p. 4111, 


lI 


15: 
CNR(dB) - P (dBW) - P (ави) 
r n (Se 
where CNR = carrier-to-noise ratio (dE) 
F = received power level (dBW) 
r 
F = receiver thermal noise level (dBW)? 
n 


ES. Digital Radio Link Parameters 


In digital systems the modem performance is 


usually plotted versus average bit energy, EM. to the 
b 
receiver noise spectral density, N. The probabilty of 
O 
error (often called bit error rate) will be determined 


from the E /N ratio. The transformation of the 
b o l 
calculated carrier-to-noise ratio (CNR) Fo ОЕ КУЧ 1 5 
b O 
written as [Ref. 23:p. 158): 


MES = το; οὔ 
E NO авт 10193108, 10log,)R 


where E = transmission noise bandwidth (Hz) 
ο 


R = transmission data rate (bit/sec) 
5.  Propability-of Bit Error Ul CULE em 

Current military multichannel communi cati ons 
use Pulse Code Modulation/Time Division Multiplexing 
(PCM/TDM) techniques to transmit digital information over 
both frequency modulated (FM) and phase modulated (FM) 
troposcatter carrier systems. The quality of the 
multiplexed circuits 15 determined by the number of bit 
errors that occur because of the channel fading and 
multipath dispersion. The probability of bit errors can 
be predicted for different FCM carrier modulation 
methods. Multiphase signaling and M-ary orthogonal 
Signaling over a Rayleigh fading channel are derived by 
J. D. Proakis [Ref. 14:pp. 490-499]. 

The bit error rate (BER) for QPSK (four-phase phase 
shift keying) and DFSK (differential phase shift keying) 15 


expressed as: 


АЕ Ер 1 2 К 
1 ти 
В πα η ο ο... у — "EM (3-14) 
Б 2 ү 2-6 k 4 E ae 
Em А 
where H = correlation coefficient 
п = (for coherent PSE) 





í sr emas (for DFS5K) 
oy 
C 
where y. = average received E /N per channel 
b o 
К maer age received E /N per bit 
b o 
L Ус 
y, = سس‎ 
J 
where L = order of diversity 
J = 1 (for BFSK signaling?) 


< (for OFSK signaling) 


1 
Eit error probabilities are depicted in Figure 3-10 for 
two-phase and four-phase DFSE signaling with L = 1,2 and 4. 
Orthogonal signaling may be viewed as M-ary 
FSK (Frequency Shift Keying). The expression for the 


probability of symbol error (F ), derived by =Froakis, 





М 
assuming no diversity (L = 1) is: 
M-1 
c m+ 1 
< ie) m (3-15) 
io = A ο 
M ES 
: 1 + + πι 
ML i Yc 
where M = “2 (for BFSE signaling) 
W = 4 (for QFSK signaling) 


The equivalent bit error rate (BER) can be computed using: 


(3-15) 
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Figure 3-10 Frobability of Error UE T MES 


(After Froakis, Ref. 14:p. 492) 


where k = MEN 

The graphs of F kagus ΕΝ tor; = and.L 1, 
M b o 

2, 4 are shown in Figure 3-11. 

The Distortion Adaptive Receiver [DAR] is 
currently being used in the military troposcatter digital 
radio set,  AN/TRC-170. Experimental modem performance 
results are illustrated in the BER versus E /N curves, 
Flgure(s) 5-12 and 5-15. The results are CA for 


three (3) different multipath delay values. 


Sunde, [Ref. 24:pp. 144-2141 has derived a 


general expression for the maximum differential 
transmission delay. This equation is valid when the 
transmitting and receiving antenna beamwiaths are 


di frerent: 


where d path distance (meters) 


0 = scatter angle (mrad) 

g 
Om transmitter take-off angle (mrad) 
Û = receiver take-off angle (mrad) 


The time dispersion (multipath spread) relative to the 
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Figure 35-13 DAR Ferformance - Dual Diversity 


(After Zawislan, Ref. 4:p. 4) 
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DAR Performance ¬ Quad Diversity 


(After Zawislan, Ref. 4:p. 4) 


average time delay for unequal antenna beamwidths 1s: 





Š Ó 
SSI PF 8 EIB) 
2С (ОХО m/s) 
в 
where c= 3x10 m/sec 


In the case of equal transmitting and receiving antenna 


τν 
` 


beamwidths, Equation 2-176 is simplified to: 


ES 


A = —— ( —17C) 
TOG 
For a partıcular order of diversity, the probability of 
error сап be computed for the AN/TRC-170 troposcatter 
radiolink Dy evaluating a -least-squares polynomial 
equation derived from experimental data curves. 
4. Fade Margin and System Reliabilty 


The system fade margin, in decibels, is the 


difference between the "practical threshold" level and 
the median received signal level. The propagation 
reliability values for the worst fading condition, 


Rayleigh fading, can be compared with their required fade 
margins in fable II ке 1B. 2-3]. The “practical 
threshold", or minimum acceptable received signal _ level, 


cannot be below the FM improvement threshold 


< 


TABLE 


RAYLEIGH FADING FROFAGATION RELIABILITY 


Single Hop 

Fropagation Reliability Fade Margin 
(7.) (dB) 
9O. O 8 
ση 18 
ταν E 28 
πο ο S 
τοσο 58 
οσο 48 





[Ref. 25:p. 711. The fade margin can be evaluated by: 


Fade (dB) = P (dBm) _ L 80) (3-18) 


where N = "practical noise threshold" 
p : 
N = Р (dBm) + FM improvement 
р n 
F = -174 dBm + 10 109 В + NF 
n 19 ΙΕ 


WU. T. Barnett and ñ. Viqants of Bell Tele- 
phone Laboratories, LRef 23:pp. 59-60], have developed an 
empirical method to determine the nondiversity annual 
path availabilty. Hacnett's procedure begins by defining 
LJ as the nondiversity annual outage probability and rc 


ndp 
as the fade occurrence factor: 


actual fade probability 
r = Rayleigh fade probabilty 
If F is the fade margin in decibels: 
actual fade probabilty 
r = 19 
10 


For the worst month: 


£ 
-5 3 
г = а х 10 --) а (5-19) 


where d = path length (statue miles) 


aA 


f = frequency (GHz) 


F = fade margin (dB) 

a = 4 (for smooth earth, over water, flat desert) 
a = 1 (for average terrain with some roughness) 

а = 9.235 (for mountainous terrain) 


Considering the annual fade occurrance: 


Г = br (3290 
yr m 
where b = б.5 (for hot, humid coastal areas) 
b = .«9 (for normal, temperate or subarctic) 
b = O.125 (for very dry climate) 


Finally the nondiversity annual path outage is: 


B = 107Ё/190 ` 


The annual nondiversity availability percentage is: 


А = 100(1 - ὑπ». (percent) 


(S= = 
The percentage of availability is improved 
with the use of frequency and space diversity. Figure 


2-14 is used to graphically determine the approximate 


availability improvement for various percentages of 
frequency separation as. o0. Figure(s) 3-15 (Dual 
Diversity) and 2-16 (Quad Diversity) provide a graphical 
method to determine the percent of path availability 
(percent of level exceeded) for different diversity 
combining techniques. 
J. Diversity Requirements 
The antenna spacing, in meters, required for 

effective space diversity has been experimentally derived 
for frequencies greater than 1 GHz, in the horizontal as 


[Ref. l1i:pp. 145-1442: 


А = 0.36 (D* + 1600) 1/2 (meters) 
(25a) 
and in the vertical as: 
2 ΟΠ 
3 y - 0.36 (D + ο... (meters) Ἢ" 
where D = parabolic antenna diameter (m) 


A satisfactory frequency separation, in MHz, for frequency 


diversity has been dervied as: 


D 
A, - (.44£/ d) (D^ & 225) 77 (mug) — (s-23c) 
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Figure 2-14 Approximate Interference Fading 
Distribution versus Order of Diversity 
and Frequency Separation 


(After het. 18:p. 2-46) 
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where + 
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transmitter frequency (MHz) 


scatter angle (mrad) 


d path length (km) 


6. Analog Radio Link Farameters 
The analog radiolink considered will be a FM 
transmitter using frequency division multiplexing (FDM). 
a. Receiver IF Bandwith 
For FM modulation techniques, the re- 


ceiver IF bandwidth is computed as [Ref. 18:p. 8-151: 


Bur 2( AF, τ (3-24) 


where Е 


p 
F 


m 


peak frequency deviation (Hz) 


maximum modulating frequency (Hz) 


The peak frequency deviation is the product of the 
frequency modulation index and the bandwidth of the 
modulating signal. The maximum modulating frequency is 
computed as the sum of the minimum modulating frequency, 
the voice channel bandwidth, and the frequency spacing 
between multiplexed supergroups. The calculated IF 


bandwidth can now be used to compute the receiver noise 


— 
/ <.) 


threshold, Equatirom “πι της carrier-to-noise ratio 15 
determined by Equation 3-12. 
b. Expected Channel Noise 
According to DCA System Performance 
Specifications (CRef. 18:p. 2-111], the channel noise 


standard for a troposcatter link is: 


L 
N(pWpO) = — FN (16,000) (X-25a) 
2000 
and 
N(dBa0) = 1010g,g(pWp) - 6 dB (2-25b) 
where G = path length (nautical miles) 
php = picowatts psophometrically weighted measured 


at, or referred to, a zero transmission level point. 

The term. dEa DENS to decibels of noise power above a 
reference noise power, with ап adjustment factor to 
compensate for equipment weighting. The referenced noise 
power that dBa 15 referred to is -85 dEm. To obtain 
dBaO, it is required to calculate the number of dB above 
this reference power the signal is. For flat voice 
channels, the corrected reference level is  —82 dEm and 


the expression for dBaO is: 


dBaO = 82 — SNR cr) 


The signal-to-noise ratio , SNR, must be calculated to 

compute Equation 3-26. The channel SNR may be calculated 
after the carrier-to-noise ratio (CNR) has been 
determined by Equation 2-11. The relationship between 
channel SNR and system CNR in a FM/FDM system is 


Breet. l:p. 2721: 


SNR = CNR + D. + FM. OE P. (oe) 
im im ip im 
where E = FM improvement factor assumed 20 dB) 
lin 
D = diversity improvement factor (Fig. 5—17) 
im - 
P = preemphasis improvement factor (Fig. 2-18) 
im 
i. = ~j(j + τω log MN (dE) 
T 10 


where M = number of voice channels 
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Figure 2-17 SNR Improvement from Diversity Techniques 


(After Freemann, Ref. 1:p. 207) 
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Figure 3-18 Freemphasis Improvement for FDM Channels 


(After Freemann, Ref. 1:p.224) 
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gain caused by the presence of an elevated troposhperic 
duct. The calculation must provide a result that was 
within 10 dE of the FDUCT computer program рей: Ес 
Seyeral related works were investigated for 4 possibie 
approach to the problem. Military troposcatter systems 
have introduced several constraints to the height gain 
problem. These constraints include: 
1. The transmitter and receiver terminals are both 
located outzide and below the elevated duct. 
2. The operating frequency range was between 4.5 GHz 
CO ο aes 
5. The geographical area of System deployment was 
limited to Western Europe (i.e. West Germany). 

A statistical analysis of elevated duct Aistograms 
was obtained from several radiosonde recording stations 
located throughout West Germany ΓΓΕΤ, 51g τοι 
information is condensed in Table III. Elevated ducts 
that had a percentage of occurrance greater than twenty 


tactically practical method of calculating the 


IV. HEIGHT GAIN COMFUTATION 


— — — Sete —  — وس جیب‎  — — —— —  —  — ων 


GENERAL 


Όπε of the purposes of this research was to develop 


height 


percent (207%) over the 3 year recording period were prime 
candidates for further study. Eight (8) elevated ducts 
were selected as typical for the area. Their optimum 


coupling heights ranged from 951 to 1452 meters above the 


surface. The elevated duct intensities were all less 
than 6 M-units with 4 M-units being the dominant value. 
TABLE 111 


ELEVATED DUCT HISTORICAL INFORMATION 


Station Mean Optimum Mean Duct Mean Duct 
Location Coupling Height Intensity Thickness 
(LAT/LONG) (Meters) (M-Units) (Meters) 
Stuttgart, FRG 

(48-49N/69-11E) 1455 το Bro 
Essen, FRG 

(51-23N/06-58E) 1412 3 106 
Hannover, FRG 

(S2-2EN/O9-41E) 1575 4 τος 
Rheine, FRG 

Ec 15M/07-25E) 1237 € 115 
Idar-Überstein, FRG 

σι ο 15) TISL š 83 
Emden, FRG 

(S3-22N/07-13E) 1071 = 124 
Goch, FRG 

(31-460H/O6-10E) 1044 & Iu 
Greifswald, FRG 

ESd—05N/1-—27ZE) Dt 4 tis 


(After 


Ortenburger, 


Ref. “£6, 


Height gain values, in decibels, were computed by the 


FDUCT program for each of the eight selected ducts. The 
following common input parameters were entered into the 
FDUCT program for each duct. ГКе+. ЖШ 


Ds The transmitter site elevation was fixed as the 


= 


reference at 5 meters above the surface. 


3 


a The receiver site elevation was inecreatea TET 


c 


J. to 280 meters in Z5 meter increments. 


The path distance was increased from 75 to 325 
kilometers in 5ο meter increments. 


D 
— 


t for both БЕКЕ 


5 пе antenna polarization waz se 
zontal and vertical polarization. 

о. fhe frequency of interest Was increased troi 
4500 MHz to 5000 МНЕ іп 100 Mhz inert enenes ee 
frequency was changed every FDUCT program run 
for each selected duct of interest. 

& The relative permittivity (15), conductivity 


(9.91 mho/m), and maximum mode attenuatian 
{1.0 dB/km) were selected for the path character— 
istics. 


B. MATHEMATICAL FORMULATION 


The FDUCT height gain results for each duct were 


investigated, and the following observations were 
discovered. The height gain values increased 
exponentially as the receiver elevation approached the 


elevated duct. The other 


ul 


Gttom OF т 


[T 


Ii 


ignificams 


f 
ГЕ 


observation was th for any fixed receiver elevation, 


the height gain increased linearly with increased path 


distance. 

Two empirical methods were considered in developing 
the height gain prediction model. One method was to 
directly το the individual height gain results for 
each of the elevated ducts and then to interpolate а 
ferent gdin result from a database. This would have 
required an extensive height gain database. The chosen 


approach was to approximate height gain curves from the 


FDUCT results using a least-squares curve їл ЕПС 
program. A second order polynomial equation was derived 
for each frequency of interest at the initial 25 


kilometer path length for each selected duct height. 

The height gain program module consisted of eight (5) 
optimum coupling height decision regions which 
E сол еа to the zelected elevated ductsz. Each duct 
height region contained a height gain polynomial equation 
for each frequency of interest. Far a particular 
Operating frequency and optimum coupling height, а 
baseline height gain was estimated and multiplied bv an 
incremental range correction factor. This range Factor 


VJ e 


ul 


the average differential difference in height gain 
between the successive SO kilometer path increments. 
Figure 4-1 outlines how the height gain estimate was 
computed. 
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Figure 4-1 Herght Gain Estimatrontr Is as 


C. RADIOQSONDE DATA ANALYSIS 

Elevated ducts are described by their optimum 
coupling height, duct intensity, and duct thickness. 
This information can be determined indirectly from 
real-time radiosonde data taken along the system’s path 
length. The atmospheric pressure, vapor pressure, 
temperature and relative humidity readings are used to 
calculate the refractive index, N, for various altitudes 


using the expression, [Ref. 27:p.41: 


ZO ΗΒ 
СБ T ОЛО; 
N = — P+ — 
E ΤΙ πο MT 
where P = atmospheric pressure (millibars) 
T = temperature (celsius) 
F = saturated vapor pressure (millibars) 
Ww 
Н = relative humidity (4) 


The radiosonde information can be entered into the 
program. The refractive index values are calculated for 
available altitude readings and then converted to 
modified refractive index values. At this point a 
M-profile plot can be obtained along with a radiosonde 
data listing. If an elevated duct does exist, the 


optimum coupling height, intensity, and thickness can be 


e) 


determined directly from the M-profile. These duct 
descriptors must satisfy the respective limits 
established by the height gain prediction model. If the 
detected elevated duct exceeds the limits of the model, 
the height gain prediction will become inaccurate, The 
program user will be alerted of “nie condition. At Chee 
Point the design engineer must decide to obtain the 
height gain from an alternative computation or neglect 


ducting effects. 
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A. GENERAL 


For validation purposes, the TROFO Frogram was used 


to determine the design predictions for a typical 
tropospheric scatter communications system. Table IV 
outlines the proposed system specifications and 
assumptions. Both the terrain and radiosonde data have 
Deen assumed to accommodate the program's capabilities- 


The program results are illustrated in^ Figure(s) ae 
through 35-4. The Radiosonde Environmental Data Listing 


(Figure 3-3) identifies the occurrence of a tropospheric 


duct by detecting a "trapped" radio ray path condition. 
Gthar refractive bending conditions identified ars 


Ssuper-refractive (bending toward the e@arth’s surface:, 
normal (standard bending), and sub-refractive ‘uoward 
Memeing). [he M—-Frofile Plot, Figure 5-4, graphically 
Verifies the program's computed duct parameters. 
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elevated duct parameters. Height gain curves were 
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TABLE IV 


FPROFOSED EXAMFLE SYSTEM SEECIFIG Bin = 


Site: Transmitter Receiver 
Latitude: Gg ο ου S7 ΟΠ J 
Longitude: UG E> ος B OS εν 
Elevation: 13 meters ^ 94 meters 


Terrain Data Available: 


Radio Terminal: Military 
Diversity: Dual 

Frequency: 4550.0 MHz 
Ete 
Ft: 


Antenna Diameter: 15 
Waveguide Length: 25 


Digital Trunk Data Rate: 


Transmission Bandwidth: 


Radiosonde 
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Data Available: 


Near Obstacle Path Mode 


(ANATRC = PONV) 


Farabolic 
per Antenna 


2.048 MbZs 
MHz 
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Gurface Duct Detected 


25 


plotted for two test cases: (1) a fixed optimum coupling 


height and range with a varying frequency, and (2) a 


ry 


fixed optimum coupling height апа frequency with 


co 


varying range. Figure(s) S-Sa thru S-Sc has shown that 


f 


prediction error of approximately 5 dB is possible. The 


error increases as the frequency approaches the adjacent 


frequency increment. In this case height gain values 
validated at 4700 MHz have been averaged inta those 
computed for 44600 MHz to produce a shift at the 
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EEedi:iction curve. Figure(s) 5—ба thru S3-éc il 


the results of onlv changing the path distance. In this 
case the greatest error detected was within 3 dB πὶ 


FDUCT, which Supports the linearity between height 
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ΟΠΕ TRANSMITTER RECEIVER 
LATITUDE: 52.00 М 5.15 N 
LONGITUDE: Gr 1 9.90 h 
ELEVATION: 15.00 M 94,00 M 


TERRAIN PROFILE TYPE: NEAR OBSTACLE PATH MODE 


TRANSMITTER TAKE-OFF ANGLE: 4.73 MRAD 
RECEIVER TAKE-CFF ANGLE: 6.97 MRAD 
SCATTER (ANGULAR DISTANCE): 20,94 MRAD 
TRANSMITTER TAKE-OFF ANGLE: cel ματ 
RECEIVER TAKE-OFF ANGLE: soo JESRESS 

SCATTER (ANGULAR DISTANCE): 1.20 DEGREES 
TRANSMIT FREQUENCY: 1560.00 N37 
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Figure 5-1 Example System Design Specifications 
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Figure 5-2 Example System Terrain Frofile Flot 
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Figure 5-3 Radiosonde Environmental Data Listing 
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Figure S-Sa Height Gain Model Comparison 


(Frequency: 4625.0 MHz) 
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Figure 5-6b Height Gain Model Comparison 


(Frequency: 4650.0 MHz) 
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Figure S-6c Height Gain Model Comparison 


(Frequency: 4575.0 MHz) 
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Figure 5-7а Height Gain Model Comparison 


(Kange: 75 kilometers) 


(2 ©) 


0.25 
ДУ Ж ЖР БЛЕЙК A οι А 


0.20 


0.15 
A e УКЕ. 


— 
= 
x 
т 
© 
L.J 
= 


0.10 
ντε. OPS 07: 


0.05 


—40.0 -55.0 
HEIGHT GAIN (DB) 
H:1.257KM F:4990MHZ R:f/5KM 


k Ë "T YES eee oque menm cna rA b E i a e D a a den ден 





Figure 5-7b Height Gain Model Comparison 


(Range: 175 Kilometers) 
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Figure 5-7c Height Gain Model Comparison 


(Range: 275 kilometers) 
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VI. SUMMARY 


A. CONCLUSIONS 

An interactive, computer-aided design program far 
military tropospheric scatter communications systems was 
developed. The program provides the communications 
engineer with an efficient and acccurate system 


performance algorithm that can be used in a tactica 


fend 


environment on A microcomputer. The following thesis 
Gb jectives were accomplished: 


1, Both analog and digital troposcatter radio 
terminals are included as design options. 


2. The influence of elevated tropospheric ducting an 
systems installed throughout West Germany was 
considered. A height gain estimation model was 
derived using numerical results obtained from the 
PDUCT main-frame computer program. 


ος System performance predictions can be calculated 
for the newly deploved digital troposcatter radio, 
AN/TRO-170, 
4. Real-time radiosonde information can be analyzed 
to determine the presence of both surface and 
Glevated tropospheric ducts. 
The accuracy of the height gain model was bounded Dy 
the range of selected duct parameters. However height 


gain results of elevated ducts which satisfied the limits 


of the model were within 3 dB of the FDUCI conus 


model. The model did not compute the height gain for 


zur face-based ducts. 


E. RECOMMENDATIONS 
Addtional research is required and should focus 
on the following areas. 


1. The validation of the computer-derived system 
perfomance results can be conducted during 
standard operational testing. Real-time radiosonde 
data should be obtained from available tactical 
weather facilities and tropospheric ducting 
phenomena determined. Received signal strength 
data can be recorded and statistically analyzed. 


2. The height gain model should be expanded to 
included unlimited duct parameters, e.g. surface- 
based ducts and elevated surface ducts. 


5. The effects of tropospheric ducting can be studied 
further within the laboratory. Radiated eneray, 
produced by laser light, can be transmitted 
through a fluid medium containing light-scattering 
particles. A common scatter volume can be formed 
by mirror-like apertures at the laser sources. 

The refractivity index of the medium could be 
controlled to Simulate ducting conditions. The 
received energy could be detected and experimental 
results studied. 


APFENDIX A 


The following prediction formulas were developed by 
the National Bureau af Standards (NBS) for calculating 


the long-term median basic transmisssion loss, L 


[Ке+. 12:р. 289] and CReft. 18:06. 856 επι». ы 
1. Long-term Median Easic Transmission Loss 
br e, = 9019610 - 2010315 + F(@d) - ВЕ T H am 
where f = frequency (MHz) 
d = mean sea level arc distance (km) 
F(@d) = attenuation function (dE) 
F = secattering-efficiency term (dB) 
"j = frequency-gain function (dE) 
ο : 
2. Attenuation Function 
For a particular symmetry factor, S. and 
approximate surface refractivity, No, the following 


= 


Figures are used to determine the attenuation  tactor. 


For a Surface refractivity: 


ων. "ΡΥ, Es - Refer to Figure ñ-i 


М = SOL апааа... Refer to Figure A-2 
= 

ΠΥ...» x Refer to Figure A~s 
= 

ο a Sw. 9s Refer to Figure A-4 


Scattering-Efficency Term 


The following equation can be used: 


Ns 





Е. _ = = (Ans! 
EC 1.0686 ( ΤᾺ ) (ho Ay һу, hy) 
where h = height of transmitter/receiver antenna beam 
O 
intersection (km) 
H = height from obstacle elevation baseline ta 
1 
the antenna beam intersection (km) 
h = transmitter obstacle elevation (πὶ 
it 
h = receiver obstacle elevation (m) 
Ir 
4. Frequency-Gain Function 
The frequency-gain function is expressed by 
as: 
το. um LH κου 
"e OE - H (ap) =) 
О О 
2 
and r = 41.92 0 tf h 
Í te 
r = 41.92 O£ h 
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for N 


Function FC @d) 


Figure Α-4 


(After Ref. 
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te 
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scatter angle (angular distance) (mrad) 
frequency (MHz) 
transmitter elevation (km) 
receiver elevation (km) 
а / B (symmetry factor) 
О O 
r /(s)(r ) 


f 1 
(r ) are graphically derived using 
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